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ABSTRACT 

We present a detailed analysis of the complex absorption apparent in the 2-6 keV X- 
ray spectrum of the bright nearby Seyfert galaxy NGC 4151. We first utilize the large 
bandpass and medium spectral resolution afforded by BeppoSAX data to construct a 
1-100 keV spectral template, which assumes the absorption arises in both warm (i.e. 
partially photoionized) and cold gas present in the line of sight to the active nucleus 
of the source. Application of this spectral model to an ASCA "long-look" observation 
of NGC 4151 reveals a partial correlation between the underlying continuum flux and 
the ionization state of the warm absorber. Such a correlation is an intrinsic property of 
a warm absorber and argues strongly in favour of this interpretation for the complex 
absorbing column over alternative partial covering models. 

The warm absorber in NGC 4151 has a column density of ~ 2 x 10 23 cm~ 2 
with ionization parameter in the range log(£) ~ 2.4 — 2.7. The inferred relatively 
low density (~ 10 5 cm~ 3 ) for the warm gas, implies an equilibration timescale for 
the dominant ions of the same order or longer than the timescale of the continuum 
variability. It follows that the warm component will invariably be observed in a non- 
equilibrium ionization state. The warm absorber in NGC 4151 may be identified as a 
multi-temperature wind produced by evaporation from the inner edge of an obscuring 
torus as discussed in a recent paper by Krolik & Kriss (2001). The unusually complex 
character of the absorption seen in NGC 4151 may then be explained in terms of 
a fortuitous line of sight which grazes the top edge of the obscuring torus so as to 
intercept a substantial column of both the warm and cold gas. 

We also find that (i) the reported hardening of the spectrum of NGC 4151 as 
the continuum level falls may be simply due to the presence of an underlying (hard 
and relatively constant) Compton-reflection component and (ii) the iron Ka line has 
a relatively narrow Gaussian profile and a line flux that remains constant over both 
short (days) and long (months to years) timescales - a relativistically broadened iron 
Ka feature was not required in our modelling. 

Key words: galaxies: active - galaxies: Seyfert - X-rays: galaxies - galaxies: NGC 
4151. 



1 INTRODUCTION 



The Seyfert 1 galaxy NGC 4151 was identified as an X-ray 
source over thirty years ago (Gursky et al. 1971). As one 
of the brightest Active Galactic Nuclei (AGN) accessible in 
the X-ray band, it has since been extensively studied by mis- 
sions such as EXOSAT, Gmga, ROSAT, ASCA, CGRO and 
more recently RXTE, BeppoSAX and Chandra. This obser- 
vational focus has revealed that the spectrum of NGC 4151 
from 0.1-100 keV is comprised of a complex mixture of emis- 
sion and absorption components, probably originating in a 
variety of locations from the innermost parts of the putative 
accretion disk in NGC 4151, out to the extended narrow- 
line region of the galaxy. Given this spectral complexity it 
remains to be seen whether NGC 4151 should or should 



not be considered as an archetype of its class (Ulrich 2000; 
Zdziarski et al. 2001). 

The intrinsic X-ray to 7-ray continuum emanating from 
the active nucleus of NGC 4151 appears to be produced 
by the thermal Comptonization of soft seed photons (e.g. 
Haardt & Maraschi 1991; Zdziarski et al. 1994, 1996, 2000; 
Petrucci et al. 2000). There is also a contribution from re- 
processing in the form of Compton-reflection and iron-K flu- 
orescence components (e.g. Maisack et al. 1991; Yaqoob et 
al. 1995; Zdziarski et al. 1996; Warwick et al. 1996; Piro et 
al. 2002; Zdziarski et al. 2001). Below ~ 5 keV the hard 
continuum is strongly cut-off by photoelectric absorption in 
a substantial (Nh ~ 10 23 cm -2 ) line-of-sight gas column 
density (e.g. Holt et al. 1980; Yaqoob et al. 1993; Weaver 
et al. 1994a,b). One of the long-standing problems in X- 
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ray studies of NGC 4151 has been to understand the ex- 
act nature and origin of this absorption. X-ray observations 
show very clearly that the cut-off is not as abrupt as one 
might expected for absorption in a uniform slab of cold so- 
lar abundance material (i.e., the increase in opacity with 
decreasing photon energy is less rapid than predicted for 
a uniform cold absorber). Various solutions have been sug- 
gested including an inhomogeneous cold absorber (the par- 
tial covering model; Holt et al. 1980; Weaver et al. 1994b), 
an absorber with grossly non-solar abundances (Yaqoob & 
Warwick 1991) and the partial photoionization of the ab- 
sorbing medium (i.e., the warm absorber model, e.g., Kro- 
lik & Kallman 1984). However, none of these has provided 
a really satisfactory explanation of the nature of the ab- 
sorber or the processes which give rise to the large changes 
in the absorption apparent in NGC 4151 on timescales of 
days or longer (Yaqoob et al. 1989, Yaqoob et al. 1993). Un- 
fortunately the presence of additional soft X-ray emission 
components which first appear at ~ 2 keV and dominate 
the spectrum below ~ 1 keV (Weaver et al. 1994a; War- 
wick, Smith & Done 1995) adds further to the complication. 
EINSTEIN and ROSAT HRI measurements (Elvis, Briel 
& Henry 1983; Morse et al. 1995) and more recently ob- 
servations by Chandra (Ogle et al. 2000; Yang, Wilson & 
Ferruit 2001) have revealed that much of this soft emission 
emanates from a spatial resolved (~ 1.6 kpc) highly ionized 
plasma coincident with the optical narrow-line region of the 
galaxy. 

In the present paper we investigate the nature of the 
complex X-ray absorption column in NGC 4151 and in par- 
ticular focus on the question of whether at least some of the 
properties of the absorber can be explained in terms of pho- 
toionization effects. Our approach has been to use a high 
signal-to-noise BeppoSAX observation to define the param- 
eters of a "spectral template" representative of the under- 
lying continuum and other features in the X-ray spectrum 
of NGC 4151. We then use this template to investigate the 
spectral variations which occurred in NGC 4151 during a 
"long-look observation" carried out in the latter stages of 
the ASCA mission. The plan of the paper is as follows. In §2 
we specify our spectral template model and determine the 
relevant spectral parameters using the BeppoSAX measure- 
ments. In §3 we give brief details of the ASCA long-look 
observation, outline the data reduction techniques employed 
and characterize the flux and spectral variability apparent 
in the source light curve. In the next section we describe the 
detailed spectral modelling of the ASCA data focusing on 
the evidence for a photoionized absorber. Finally in §5 we 
discuss the implications of our results. 



2 THE BeppoSAX DATA AND THE SPECTRAL 
TEMPLATE 

The broad bandpass of the instruments on BeppoSAX is a 
tremendous help when attempting to constrain the form of 
the hard X-ray spectrum of NGC 4151. During the period 
1996-1999, this source was observed three times with Bep- 
poSAX as reported by Piro et al. (2002). Here we utilize the 
observation carried out in January, 1999, which provides the 
best signal to noise ratio of the available datasets. X-ray data 
from three of the four instruments on BeppoSAX (namely 



the LECS, MECS and PDS instruments - see Parmar et 
al. 1997, Boella et al. 1997, Frontera et al. 1997 respectively 
were obtained from the ASI Science Data Center. In the case 
of the LECS and MECS instruments, source spectra were 
preprocessed via standard BeppoSAX procedures using on- 
source extraction regions of 6' and 4' radius respectively. 
Standard blank-sky background files were available for the 
LECS and MECS whereas in the case of the PDS instru- 
ment, background-subtracted source spectra were supplied 
directly. The spectra were accumulated over the full obser- 
vation interval (see Piro et al. 2002) to give an on-source 
exposure time of 83 ks in the MECS and roughly half this 
for the other instruments. 

A BL Lac object (MS 1207.9+3945) is located ~ 5' to 
the north of NGC 4151. Piro et al. (2002) note that its count 
rate is ~ 2% and ~ 20% of that of NGC 4151 in the MECS 
and LECS instruments respectively. The contamination of 
the NGC 4151 spectra will be negligible for the MECS but 
may, on the basis of the Piro et al. (2002) figures, be as 
high as ~20% in the LECS at ~ 1 keV. In practice the 
presence of this confusing source in the LECS data will have 
little overall impact on our spectra fitting analysis (given the 
relatively soft spectrum of the BLLAC - Warwick, Smith & 
Done 1995). 

Spectral fitting was carried out using the XSPEC 
Vll.0.1 software package with the Sept 1997 instrument 
response matrices. For this purpose the source spectra were 
grouped to give a minimum of 20 counts per bin and the 
data below 1 keV (from the LECS) were excluded to avoid 
complications relating to the form of the soft X-ray emitting 
components and presence of the soft confusing source. We 
adopt a spectral template model which includes the follow- 
ing emission components: 

(i) A power- law continuum with a normalization, Ai, and 
photon index, F, exhibiting a high-energy break at 100 keV; 

(ii) A neutral Compton-reflection component (modelled 
by PEXRAV in XSPEC, Magdziarz & Zdziarski 1995) with 
only the reflection scaling factor, R, as a free parameter. The 
parameters relating to the incident continuum were tied to 
those of the hard power-law component. In addition cos i 
was fixed at 0.5 and the metal abundance in the reflector 
was fixed at the solar value; 

(iii) An iron Ka emission line of intensity Iku at an en- 
ergy Ekc (with an intrinsic line width ok a set to 0.1 keV); 

(iv) A second power-law continuum representative of all 
the additional soft X-ray emitting components with a slope 
coupled to that of the primary hard continuum but with a 
free normalization, Ai. 

As noted earlier the complex absorption in NGC 4151 
has in the past been modelled in number of ways, with ar- 
guably the most successful approach being the partial cov- 
ering (sometimes referred to in the literature as a dual ab- 
sorber) model (e.g. Holt et al. 1980; Yaqoob & Warwick 
1991; Weaver et al. 1994b). In the partial covering scenario a 
fraction f cov of the hard continuum is absorbed by a cold gas 
column density Nh,i, whilst the remaining 1 — f cov fraction 
of the continuum intercepts a reduced column density Nh,2- 
In the present paper we prefer to focus on a model employing 
the same number of free parameters but with the absorber 
stratified along the line of sight rather than perpendicular 
to it. In our preferred scenario the complex absorber is rep- 
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Table 1. The results of fitting the "spectral template" 
BeppoSAX observation 



to the 



Model 
Parameter 



Best-Fit 
Value 



Units 



^ H ,warm 

r 

R 

Ai 
M 
Eku 

x 2 

d-o-f 



23.9 
3.4 
1.65 
0.37 
2.48 
4.6 
1 

6.37 
3.4 



+0.9 
0.9 
+1.1 
0.8 
+0.02 
0.03 
+0.18 
0.15 
0.0(1 
0.03 
+0.2 
0.2 
+0.1 
-0.1 
+0.03 
0.02 
+0.4 
0.3 



10 22 cm 
10 22 cm 



10 _2 photon keV -1 
10~ 3 photon keV -1 
keV 

10 — 4 photon cm -2 s 



832 
815 



resented as product of two absorption components, namely 
a warm column density NH.warm and a cold gas column 
Nh.coU- For the former we use multiplicative table models 
in XSPEC generated via the photoionization code XSTAR 
(Kallman 2001). The ionization state of the warm gas is gov- 
erned by the ionization parameter £ = Li 0n /n r 2 where Li 0n 
is the source luminosity in the 0.0136-13.6 keV bandpass in 
erg s _1 , n is the number of hydrogen atoms/ions in the gas 
per cm -3 and r is the distance from the central source to 
the inner edge of the warm cloud in cm. We assumed an 
ionizing continuum of the form adopted by Krolik & Kriss 
(1995). For further details of the photoionization modelling 
see Griffiths et al. (1998). The cold absorption was repre- 
sented by the wabs model in XSPEC which is based on the 
absorption cross-sections tabulated in Morrison & McCam- 
mon (1983). Here we assume solar abundances in both the 
warm and cold absorber. 

The adopted spectral model assumes that only the hard 
power-law continuum is subject to the complex absorption. 
However absorption arising in the line-of-sight column den- 
sity through our own Galaxy is applied to all four emission 
components (Nh,Go,i — 2 x 10 20 cm -2 ). 

The results of fitting this spectral template to the Bep- 
poSAX data are detailed in Table [j], which lists the best 
fitting values for the nine free parameters of the model. The 
quoted errors (here as elsewhere in this paper) are at the 
90% confidence level as defined by a Ax 2 =2.71 criterion 
(i.e. assuming one interesting parameter). 

The spectral template provides an excellent fit to the 
BeppoSAX data as illustrated in Fig. |l|. Fig. |^ shows the 
corresponding best-fitting model spectrum. 



3 THE ASCA LONG-LOOK OBSERVATION 

In the present paper we concentrate on the final observa- 
tion of NGC 4151 carried out by ASCA during the period 
May 12-25, 2000. This was an exceptionally long observa- 
tion carried out during the last phase of the ASCA pro- 
gramme, fully meriting its description as a "long- look". The 
ASCA payload included four X-ray telescopes, two equipped 
with Solid-state Imaging Spectrometers (SIS) and two with 
Gas Imaging Spectrometers (GIS). All the relevant datasets 




Figure 1. The BeppoSAX LECS (black), MECS (red) and PDS 
(green) data fitted by our spectral template model. Upper panel: 
The count rate spectra and best-fit model. Lower panel: The ratio 
of the data to the best-fit model. 




Energy (keV) 

Figure 2. The spectral template model giving the best-fit to the 
BeppoSAX data. There are four emission components comprising 
of a hard power-law continuum (red), Compton reflection (green), 
an iron Ka line (dark blue) and a second power-law represent- 
ing soft X-ray emitting components (light blue). Complex absorp- 
tion (warm+cold components) causes the low-energy cut-off in 
the hard continuum. All the emission components are subject to 
Galactic absorption. 



pertaining to the SIS and GIS instruments were obtained 
from the ASCA public archive at the HEASARC. Standard 
screening criteria and reduction techniques were employed 
using software routines within the FTOOLS package. During 
the long-look observation the SIS instruments were operated 
in a single CCD mode with the nearby BL Lac object posi- 
tioned at the very edge or just off-chip in the SIS-0 and SIS-1 
instruments respectively. In an attempt to reduce both the 
impact of the instrument background and the contamina- 
tion from the BL Lac to negligible proportions, we applied a 
1' radius source extraction circle to the SIS data. Compared 
to a more standard 3' radius extraction cell, this reduced 
the 1-10 keV SIS count rate by ~ 11%. For the GIS data we 
used a 3' radius source extraction region and a background 
spectrum taken from an off-source region in the GIS field of 
view. 

As discussed recently by Turner et al. (2001) the 
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ASCA SIS detectors have shown a degradation in efficiency 
at lower energies, which is probably due to an increased 
dark current levels and decreased charge transfer efficiency 
(CTE), producing SIS spectra which diverge from each other 
and from the GIS data. Furthermore, data from the last 
phase of ASCA operations (AO-8) have revealed a non-linear 
evolution of the SIS CTE. Here we have applied the interim 
solution released on 2001 February 13 by the ASCA GOF 
(in the form of the CTE file sisph2pi-13Q201.fits) which 
reduces some but not all of the spectral inconsistencies (see 
below). We restrict our analysis to the 1-10 keV band data 
to avoid the worst of the low energy calibration uncertain- 
ties in the SIS instruments (e.g. see Appendix A of Weaver 
et al. 2001). For subsequent analysis, the total accumulated 
on-source exposure was 285 ks in the SIS and 435 ks in the 
GIS detectors. 

The light curve in the full 1-10 keV band is shown in 
Fig. ^ for the combined SIS instruments. An almost iden- 
tical light curve is obtained for the combined GIS instru- 
ments demonstrating that, at least in terms of the broad 
band variability, the two SIS and GIS instruments produce 
very similar results. During the ~ 13 day observation the 
source shows a significant brightening by a factor of more 
than 3 followed by a comparable decline. The fastest vari- 
ation is a flux increase of ~ 35% in roughly 12 hrs. Fig. ^ 
also shows a plot of the 2-5 keV/5-10 keV softness ratio 
versus time (for simplicity we show only the SIS data). The 
brightening of the source is clearly associated with a sig- 
nificant spectral softening, although the spectrum remains 
relatively soft during the subsequent decline in flux. 

We have also investigated the spectral variability ex- 
hibited by NGC 4151 by calculating the fractional rms vari- 
ability amplitude as a function of energy (see Edelson et 
al. 2002). The result, shown in in Fig. ^| demonstrates, in a 
model independent fashion, the steep rise in the amplitude 
of the variability above 1 keV to a maximum at ~ 3 keV, 
followed a significant decline at higher energy. There is also 
a hint of a dip at the energy of the iron Ka line (6.2-6.6 
keV). 



4 FITTING THE SPECTRAL TEMPLATE TO 
THE ASCA DATA 

In order to study the precise nature of the spectral variabil- 
ity exhibited by the source, we have broken the long-look 
observation up into eight sections utilizing gaps in the light 
curve to define section boundaries as illustrated in the lower 
panel of Fig ^. Appropriate response matrices and ancil- 
lary response files were constructed using standard FTOOLS 
routines and the spectra were again binned to a minimum 
of 20 counts per bin to allow the use of the \ 2 statistic in 
the spectral fitting. 

Here our goal is to explain the softening of the spectrum 
as the continuum rises and the subsequent lack of spectral 
evolution during the continuum decline in terms of photoion- 
ization effects. In the context of the spectral template de- 
scription discussed earlier, this implies changes in the ion- 
ization of the warm absorber induced by variations in the 
underlying continuum flux from the central source. 

We adopt a very stringent requirement in spectrally fit- 
ting the ASCA data in that we fix all but three of the spectral 
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Figure 3. Upper panel: the combined SIS 1-10 keV X-ray light 
curve from the "long-look' ASCA observation. Lower panel: The 
variation the 2-5 keV/5— 10 keV softness ratio during the obser- 
vation. The horizontal bars illustrate how the observation was 
split into the eight time segments for the spectral analysis. The 
coloured horizontal bars indicate the segments for which spectra 
are shown in Fig. ^. 




Energy (keV) 

Figure 4. The fractional rms variability amplitude as a function 
of energy calculated for the entire ASCA long-look observation 
(SIS data). 



parameters of the template model at the values derived from 
the BeppoSAX observation (Table |l]). The three parame- 
ters permitted to vary are the normalization of the hard 
power- law continuum, A\, the normalization of the soft X- 
ray power-law, A%, and the ionization parameter, £. The 
column densities of the warm and cold components are kept 
fixed, as is the slope of the underlying continuum and the 
parameters of the iron Ka line. For the Compton reflection, 
we fixed the actual component flux at the BeppoSAX level 
(rather than fixing the scale parameter R). 

Before proceeding with the ASCA spectral analysis we 
need to investigated the degree of inconsistency between the 
SIS and GIS instruments. A comparison of the residuals 
when the above "template" model is fitted to the SIS-0/SIS- 
1 and GIS-2 spectral data from time segment 3 is shown in 
Fig. |5| As noted earlier these correspond to data extracted 
with a 1' radius source cell in the case of the SIS and a 
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Figure 5. A comparison of the best-fitting spectral template 
model to the combined SIS (black) and G1S-2 (red) spectra from 
time segment 3 of the long-look observation. Upper panel: The 
count rate spectra and best-fit model. Lower panel: The ratio of 
the data to the best-fit model. Note the discrepancies between the 
SIS and GIS are most evident below ~ f .5 keV and above 8 keV. 



3' radius cell for the GIS. Discrepancies between the two 
datasets are most evident below 2 keV and above 8 keV. 
The contamination of the GIS spectrum by the BL Lac ap- 
pears to account for most of the problem at low energies. In 
fact we find that GIS data extracted within a 1' radius cell 
(thereby greatly reducing the BL Lac contribution) actually 
agree reasonably well with the SIS data, although this leads 
to a ~ 60% loss of counts in the GIS. The problem at the 
high energy end is similar in character to that reported by 
Turner et al. (2001) in a comparable ASCA long-look ob- 
servation of the narrow-line Seyfert 1 galaxy Akn 564. Here 
we take a somewhat different approach to that employed by 
Turner et al. (2001). In the subsequent analysis we make 
use solely of the spectra derived from the SIS data. For our 
investigation the exclusion of the soft BL Lac flux is a useful 
advantage which outweighs the spectral uncertainty above 
8 keV, since we employ a very tightly constrained model of 
the hard continuum (via the adopted spectral template). In 
any case, our main interest is to characterize the spectral 
changes rather than to determine precise values for spectral 
parameters. Consistency checks show that very similar re- 
sults to those reported below for the SIS are obtained with 
the GIS instruments (in fact somewhat better reduced \ 2 
are recorded) provided an extra soft component representa- 
tive of the BL Lac contribution is including in the spectral 
modelling. 

Fig. H compares the SIS spectra from time segments 1, 
3 and 5 of the long-look observation during which the source 
count rate in the 1-10 keV band increased by over a factor 
of 3. These SIS spectra illustrate that the level of the hard 
continuum rises by a factor of ~ 2 and that it is the marked 
softening of the spectrum in the 2-6 keV band that accounts 
for the additional factor in the count rate increase (c/. Fig. 

S- 

The results of fitting the constrained template model 
to the spectra from the eight segments of the long-look ob- 
servation are given in Table ^. Overall the model provides 
a reasonable fit to the eight spectral datasets with a com- 
bined x 2 of 4616 for 3665 d-o-f (de grees of freedom). Inter- 
estingly the spectrally variation revealed by the long-look 
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Figure 6. The 1-10 keV combined SIS spectrum of NGC 4151 
as measured in three time intervals during the ASCA long-look 
observation. 
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Figure 7. The 1—10 keV X-ray spectrum of NGC 4151 as mea- 
sured by ASCA in the period 1993-1995. These are count rate 
spectra from the SIS-0 detector only. 



data largely encompasses the range of spectral form exhib- 
ited by NGC 4151 in earlier ASCA observations (Fig. |). 
In fact our spectral template model provides quite a good 
description of all the available datasets from the previous 
seven years of observation by both ASCA and BeppoSAX , 
with a range of ionization parameter not dissimilar to that 
derived here. 

If we relax some of the parameter constraints then fur- 
ther improvements in the fit to the ASCA long-look seg- 
ments can be obtained. For example, if we allow Ekc and 
Ikc to vary (but with the parameter value tied across the 
eight spectral datasets), then we obtain X 2 = 4552 for 3663 
d-o-f, with best fit values of E Ka = 6.395 ± 0.012 keV and 
Ikc = 2.2^Q 3 x 10 _4 photon cm -2 s _1 (which according to 
the F-test is significant at greater than 5 sigma). As a fur- 
ther step if we allow the two column density components to 
vary (with respect to the BeppoSAX values) then \ 2 = 4295 



N H ,coid = 4.6 ± 0.2 x 10 22 cm -2 (which is very highly sig- 
nificant according to the F-test). In all these fits the largest 
residuals are in range 1-2 keV, where the simple power-law 
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Table 2. Spectral fitting of the eight time segments from the 
ASCA long-look observation. 



Segment log(g) 



9 zlCQ+0-009 
z.<ioo_ oog 

n c-io+0.007 
z.dlo_ 005 

9 kc: 4 +0.006 

9 61 5+0.004 
Z.01D_q q 13 

2 644+ - 001 

2 636+ 003 
Z - OJD -0 .009 

9 630+ ' 004 

o 643+0-006 

z,u,J -0.008 



1 

2.58 
3.79 
4.59 
5.25 
3.50 
2.42 
2.51 



+ 0.04 

-0.04 
+0.04 
0.04 
+0.03 
0.05 
+0.10 
0.01 
+0.08 
0.01 
+0.05 
0.03 
+0.03 
0.02 
+0.04 
0.03 



A 2 



n on+0.03 
U - 8 ^-0.03 
-,+0.02 
-0.02 
n ocr+0.02 
U - 85 -0.03 
n ocr+0.07 
U -° a -0.02 
7+0.06 
-0.01 

n s«+°- 03 

qs+0 03 
u - 9J -0.03 

94+ 04 
U,M -0.04 



0.80 4 



0.87 1 



10" 2 photon keV" 
10~ 3 photon keV" 



s" 1 at 1 keV. 
s _1 at 1 keV. 



model of the soft emission components is clearly a gross 
over-simplification (Ogle et al. 2000). However, the presence 
of these 1-2 keV residuals is very unlikely to influence our 
interpretation of the observed spectral variations in terms 
of a variable warm absorber (see below). 

The spectral fitting of the ASCA long-look observation 
reveals a range of ionization states for the warm absorber in 
NGC 4151 (Table ^), A clear correlation is apparent between 
the ionization parameter (£) and the level of the hard con- 
tinuum (Ai) for the first five segments of the observation, 
which correspond to the period of increasing flux in the light 
curve (Fig |). However as the flux drops (segments 6-8) the 
correlation is lost and the ionization state remains at a com- 
parable level to that seen at the peak of the lightcurve. In 
contrast, in the BeppoSAX observation, although the con- 
tinuum was at a level comparable to that at the peak of 
ASCA light curve, in this case the ionization of the warm 
absorber was at the bottom end of the range observed with 
ASCA (Fig||). The observed correlation in the ASCA long- 
look data lends support to the warm absorber interpretation 
for the observed absorbing column and we argue in §5 that 
the lack of a sustained correlation between ionization level 
and the continuum flux, is not necessarily inconsistent with 
this description. 

The normalization of the soft power-law component, 
A2, varies by ~ 15% during the ASCA long- look observa- 
tion. There is also an apparent factor 2 decline in the soft 
component over the 18 month interval between the Bep- 
poSAX and ASCA observations, although as noted earlier 
there may be some preferential contamination of the Bep- 
poSAX data by the nearby BL Lac object. Bearing in mind 
that much of the soft emission in NGC 4151 originates in the 
spatially extended component recently imaged by Chandra, 
(Ogle et al. 2000; Yang et al. 2001), these results would sug- 
gest that there may be some leakage of the hard continuum 
flux through the absorber down to ~ 1 keV (as suggested by 
Fig ^) whereas in the current spectral modelling the cut-off 
is nearer to 2 keV (see Fig ^ . Since the complex absorption 
in NGC 4151 is likely to arise from gas exhibiting a range 
of ionization parameters, our two component (cold + warm 
gas) model must, at best, represent only a rough approxi- 
mation to the real system. 



PowerLow Norm (x10 2 Photons keV 1 cm 2 s ') 

Figure 8. Correlation between ionization parameter (£) and the 
normalization of the hard power-law continuum for the eight seg- 
ments of the ASCA "long-look" observation. The red points mark 
the segments where the lightcurve is rising, the green points mark 
the segments where the lightcurve is falling. The point marked by 
the star symbol is from the BeppoSAX observation; the error bar 
is larger in this case due to the additional free parameters in the 
spectral fitting. 



5 DISCUSSION 

5.1 The physical state of the photoionized 
absorber in NGC 4151 

One of the major advantages of the warm absorber descrip- 
tion is that it is based on a physical model which is, in prin- 
ciple, readily testable. For example, given sufficient spectral 
resolution, sensitivity and dynamic range, the true ioniza- 
tion state of the warm gas can be determined directly from 
the absorption edges present in the X-ray spectrum, A fur- 
ther key diagnostic, somewhat more in tune with current 
measurement capabilities in the case of heavily absorbed 
sources such as NGC 4151, is that ionization parameter of 
the warm gas should track the intensity of the ionizing radi- 
ation field, albeit in a complex fashion (Nicastro et al. 1999). 
In contrast the popular partial covering (or dual absorber) 
models are based on rather arbitrary, constructed geometries 
which are very difficult to verify. 

Previous studies of NGC 4151 have failed to demon- 
strate in an unambiguous fashion that the absorbing 
medium in NGC 4151 does respond to changes in the level 
of the incident hard continuum (Yaqoob et al. 1989; Fiore et 
al. 1990; Warwick et al. 1996; Piro et al. 2002). However, the 
ASCA long- look observation provides an excellent oppor- 
tunity to study the response of the absorber to continuum 
changes. The correlation between ionization parameter and 
the level of the hard continuum apparent over the rising part 
of the light curve (Fig. ^ would appear to be the signature of 
a partially photoionized medium. However, a problem then 
presents itself in that the relatively high ionization state 
is maintained over the latter part of the observation despite 
the decline in the continuum. A possible explanation for this 
behaviour is that the timescale on which significant contin- 
uum variations are occurring (~ 3 x 10 5 s) is somewhat 
shorter than the timescale on which the balance between 
photoionization and recombination can be re-established in 
the plasma. The result is that with quasi-continuous contin- 
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uum variations (of which the light curve in Fig. g represents 
only a brief snapshot) the warm absorber is invariably in a 
non-equilibrium ionization state. Nicastro et al. (1999) have 
discussed some of the complications which arise in such cir- 
cumstances and show that the delayed response of a medium 
to sharp increases and decreases in the incident ionizing flux 
can give rise to variations in the ionization state apparently 
unconnected with flux changes (and even anti-correlations 
of ionization level with flux). In the case of NGC 4151 there 
is clear evidence that the medium is slow to respond; for ex- 
ample even during the period when the ionization parameter 
and the continuum show a correlation, the former changes 
by a factor 1.8 compared to a factor of 2.5 for the latter. 
However, the fact that the degree of ionization tracks the 
light curve better during the rising segment than during the 
decline is consistent with the fact that Ti an tx L~^ n , as is 
evident from Fig. 1 of Nicastro et al. (1999). 

The parameter values for the warm absorber in NGC 
4151 are not too far removed from one of the standard cases 
considered by Krolik & Kriss (2001) (although the assumed 
form of underlying continuum is somewhat modified with 
respect to the version in Krolik & Kriss et al. (1995), which 
we have adopted in the present work). For example, taking 
£ = 300 and Li 0n — 10 44 erg s _1 as representative values for 
the warm absorber in NGC 4151 and taking the distance of 
the warm medium from the continuum source to be 1 pc, 
then the inferred gas density is n ~ 4 x 10 4 cm -3 . Table 1 
in Krolik & Kriss (2001) gives the corresponding equilibra- 
tion timescales for selected ions in such a plasma. The warm 
absorber cut-off is best measured between 2-6 keV, where 
absorption edges due to the hydrogen- and helium-like ions 
of Mg, Si, S, Ar and Ca will be of importance (although 
not individually resolvable in the ASCA data). For these 
ions the relevant equilibration timescales are typically be- 
tween 1 — 8 x 10 6 s. In fact, the measured column density 
for the warm absorber in NGC 4151 of ~ 2 x 10 23 cm" 2 
implies r max ~0.5 pc (Krolik & Kriss 2001), which in turn 
(since ti 0n oc r 2 ), requires a factor 4 reduction in the equi- 
libration timescales quoted by Krolik & Kriss. Nevertheless 
it is still true that the variability timescale in NGC 4151 
is short enough for the warm absorber to be out of ioniza- 
tion equilibrium as required. Since different ions reach their 
equilibrium abundance, for a given flux of ionizing radiation, 
on different timescales, the detailed form of the absorption 
cut-off will depended on the luminosity history of the source 
stretching back over a period equal to the longest relevant 
equilibration timescale, which in NGC 4151 will be ~ 10 7 s. 

5.2 The properties of the iron Ka line 

In the present paper we have modelled the iron Ka emis- 
sion in a very simple fashion, namely as a single narrow 
line at an energy appropriate to cold (i.e. neutral or low- 
ionization) material. In the long-look ASCA observation, the 
measured line flux was 2.2l°; 2 x 10" 4 photon cm s , with 
no evidence for any variations in either the line intensity or 
in the line profile during the observation period. The Bep- 
poSAX observation gave a somewhat higher value (at the 3<r 
level) suggesting the possibility of line variations, but on the 
other hand the ASCA value is in excellent agreement with 
the iron Ka line flux measured by Ginga almost ten years 
earlier. The lack of significant variability on timescales of 



many years suggests the origin of this narrow component is 
in a medium located well outside of the immediate surround- 
ings of the nuclear source {i.e. not in a putative accretion 
disk), but probably does not exclude the region occupied 
by the warm absorber. However, the predicted range of iron 
ionization states in the warm medium is Fe XVI - FeXXI 
with the associated Ka line emission ranging in energy from 
6.4 - 6.6 keV. Since the measured line energy excludes much 
of this range, we conclude that the bulk of the line equiva- 
lent width (which ranges from 82 - 195 eV during the eight 
segments of the observation, depending on the continuum 
level) derives from less strongly photoionized material lo- 
cated further from the nucleus than the warm medium. The 
cold absorber in our model can only account for the ob- 
served iron Ka flux if it subtends a full 4ir steradians as 
viewed from the central source and has an iron abundance 
a factor 3 higher than solar. 

In contrast with previous claims suggesting a complex 
broad iron line structure in NGC 4151 (Yaqoob, et al. 1995; 
Wang et al. 1999, 2001; Yaqoob, et al. 2001), we find no 
compelling evidence to include a relativistically broadened 
line in any of our modelling (the best fit intrinsic line width 
of the "narrow" line in the ASCA data was 110 ± 20 eV but 
this may over estimate the true line width given the SIS cali- 
bration uncertainties) . The lack of any variability in the line 
profile during the observation contrasts with earlier claims 
of line profile variability in NGC 4151 on 10 4 s timescales 
(Wang et al. 2001). 

Clearly a very extended broad-line profile is easily con- 
fused with a continuum form including a non- negligible con- 
tribution from Compton reflection and modified by a com- 
plex absorption. The complex variable nature of the underly- 
ing continuum and the absorption may then easily manifest 
itself as variability in the broad-line profile. We will return 
to this issue at a later date utilizing newly acquired data 
from XMM-Newton. 

5.3 The spectral index versus flux correlation 

Previous observations of NGC 4151 with both EXOSAT and 
Ginga have revealed an apparent correlation between spec- 
tral index measured in the 2-10 keV band and the level 
of the underlying continuum (Perola et al. 1986; Yaqoob & 
Warwick 1991; Yaqoob et al. 1993). Fig. ^ shows a compi- 
lation of measurements from EXOSAT and Ginga (taken 
from Yaqoob et al. 1993) which suggests that as the source 
brightens the photon index steepens, from T ~ 1.35 to 
T ~ 1.7, at which point the correlation saturates. Yaqoob 
et al. (1993) also note that the spectral slope changes oc- 
cur on roughly the same timescale as the flux variations. 
However, a perplexing aspect of the spectral index flux cor- 
relation observed in the medium energy bandpass is that the 
effect is not mirrored in the hard X-ray/soft gamma regime 
(Johnson et al. 1997; Piro et al. 2002). 

An underlying assumption of the spectral template 
model adopted in the present paper is that the spectral 
index of the hard power-law continuum in NGC 4151 re- 
mains constant (at a value V ~ 1.65), as the continuum flux 
changes, in apparent contradiction to the correlation noted 
above. (Testing this assumption against the ASCA obser- 
vation demonstrates that the ASCA datasets are broadly 
consistent with the assumed value of r=1.65, however the 
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Figure 9. Measurements of the spectral photon index versus the 
absorption-corrected 2-10 keV flux from EXOSAT and Ginga ob- 
servations, taken from Fig. 2 of Yaqoob et al. (1993). The solid 
curve represents the results of a simulation in which a constant 
Compton-reflection component was included along with a power- 
law continuum with a variable normalization but fixed slope (i.e., 
r = 1.65). See the text for further details. 

values are poorly constrained due to the very limited effec- 
tive bandpass, and the fact that the values of F derived are 
very dependent on the exact details of the applied spectral 
model) . 

We have investigated the possibility that the spectral 
index versus relation may actually be a misinterpretation 
of the softening of the NGC 4151 spectrum as the source 
flux increases, bearing in mind that the reported correla- 
tion was established using both proportional counter data 
with limited spectral resolution and a highly simplified spec- 
tral model (in fact just a power-law continuum model, mod- 
ified by a heavy cold absorbing column). To this end we 
have taken our spectral template model and an appropriate 
Ginga response, and simulated a set of five Ginga spectra 
using XSPEC. In this simulation, the normalization of the 
hard power-law continuum was increased in steps from 1 to 
10 x 10~ 2 photons cm -2 s _1 , but with all the other parame- 
ter values, (including the ionization parameter and the flux 
in the Compton-reflection component) frozen at the values 
obtained from the BeppoSAX spectral fitting. The simu- 
lated 3.5-20 keV Ginga spectra were then each fitted with 
a simple absorbed power-law model ( cf. Yaqoob & Warwick 
1991; Yaqoob et al. 1993) and the effective photon spectral 
index measured. The result was a spectral index versus 2-10 
keV flux correlation closely matching the reported correla- 
tion (Fig. |^) . Since the simulation did not include any of the 
complications associated with varying ionization parameter 
of the warm absorber, it follows that the observed spectral 
changes are induced solely as a result of the presence of a 
non-varying spectrally-hard Compton-reflection component. 
Clearly, in this scenario, as the direct continuum brightens, 
the relative contribution of the Compton-reflection declines 
and the overall spectrum softens until eventually the mea- 
sured spectral index saturates at the value pertaining to the 
underlying continuum. 



6 CONCLUSION 

There has been much debate as to whether NGC 4151 is sim- 
ply a very bright, nearby representative of the broad class 
of galaxies harbouring luminous Seyfert nuclei or a unique 
object (Ulrich 2000). In the X-ray band this uncertainty is 
due, in large measure, to the heavy, complex and variable ab- 
sorption which characterizes the 2-6 keV spectrum of NGC 
4151 but which is not mirrored in the great majority of well- 
studied Seyfert galaxies (Mrk 6 is arguably the best known 
example of a NGC 4151 analogue, Feldmeier et al. 1999). 

In this paper we have presented a detailed analysis of 
two archival datasets of NGC 4151, one from BeppoSAX the 
other from ASCA, which tests the hypothesis that the com- 
plex absorption is due to a combination of warm and cold gas 
distributed along the line of sight to the active Seyfert nu- 
cleus in this source. We find that this relatively simple model 
does indeed provide a good description of the two datasets 
in question and, in fact, gives a reasonable fit to all of the 
ASCA and BeppoSAX observations from the past seven 
years. Our key finding is that the warm absorber corresponds 
to a relative low density (~ 10 5 cm -3 ) gas, with a column 
density of ~ 2 x 10 23 cm -2 and with ionization parame- 
ter in the range log(£) ~ 2.4 — 2.7. Since the equilibration 
timescale for the dominant ions is of the same order or longer 
than the timescale of the continuum variability, the warm 
component is invariably observed in a non-equilibrium ion- 
ization state. Non-equilibrium conditions help explain why 
past studies based on relatively short observations spaced by 
weeks, months or years have failed to identify a consistent 
signature of the warm absorber in this source. 

In a recent paper, Krolik & Kriss (2001) discuss how 
evaporation from the inner edge of the obscuring torus in 
an AGN can give rise to an inhomogeneous photoionized 
wind in which a broad range of temperatures coexist in equi- 
librium. Krolik & Kriss (2001) then go on to suggest that 
this wind is the origin of the highly ionized, warm absorbers 
seen in over half of type 1 Seyfert galaxies. For most Seyfert 
1 galaxies the properties of the warm absorber are inferred 
from the OVII and OVIII absorption edges, whereas in NGC 
4151 the combination of warm and cold absorbers gives rise 
to a sharp spectral cutoff below ~ 2 keV, thus eliminating 
the oxygen features. Nevertheless, as noted earlier, the in- 
ferred properties of the warm component in NGC 4151 are 
not too far removed from one the standard cases detailed by 
Krolik & Kriss (2001). It is therefore reasonable to hypothe- 
sise that the warm absorber in NGC 4151 originates in such 
a multi-temperature wind, but that the unusually complex 
character of the absorption is due to a line of sight which 
grazes the top edge of the obscuring torus so as to intercept 
a substantial column of both the warm and cold components. 
In essence in terms of its X-ray absorption properties, NGC 
4151 is intermediate between the Seyfert galaxies in which 
the absorption, if present, is predominantly due to the Kro- 
lik & Kriss's multi-temperature wind (type 1 objects) and 
those where the cooler material of the obscuring torus dom- 
inates (type 2 objects). 

The other findings of this paper are that (i) the reported 
hardening of the spectrum of NGC 4151 as the continuum 
level falls may be simply due to the presence of an underly- 
ing (hard and relatively constant) Compton-reflection com- 
ponent and (ii) there is no compelling evidence for a rela- 
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tivistically broadened iron Kq line in NGC 4151 contrary to 
earlier claims. 

In summary, it appears that many of the X-ray prop- 
erties of NGC 4151 can be explained in terms of the cur- 
rent paradigm for Seyfert galaxies. Future observations with 
XMM-Newton, Chandra and INTEGRAL will be able to test 
the degree to which our spectral template description is valid 
and should further illuminate the question of whether NGC 
4151 should be regarded as the archetypal Seyfert galaxy. 
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